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Abstract 

Just as for the ordinary quantum harmonic oscillators, we expect the zero-point energy to play 
a crucial role in the correct high temperature behavior. We accordingly reformulate the theory of 
the statistical distribution function for the q-deformed boson oscillators and develop an approxi- 
mate theory incorporating the zero-point energy. We are then able to demonstrate that for small 
deformations, the theory reproduces the correct limits both for very high temperatures and for very 
low temperatures. The deformed theory thus reduces to the undeformed theory in these extreme 
cases. 
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Let us begin with the classical distribution function of ordinary quantum boson harmonic oscil- 
lators, when the zero point energy is discarded. The occupational probability 



and Z the canonical partition function are given by the well-known results 
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where x = (3hu. The distribution function is 
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the familiar result for the ordinary harmonic oscillators. The q-deformed theory is described by the 
g— analogue of the algebra of creation and annihilation operators of the boson harmonic oscillators, 
H - ] defined by 

[a, a] = [a^, a^] = 0, aa^ — qa^a = (4) 

where q is the deformation parameter and N is the number operator which obeys the commutation 
relations 

[N,a)] = a\ [N,a] = -a. (5) 
The distribution function is [H the g— generalization of Eq.@, thus 

/. = E H Pn (6) 

where [n] stands for Uq and is the set of basis numbers, defined by 

This is in contrast to the case of the ordinary boson oscillators wherein the summation involves 
n, the eigenvalue of the number operator. Accordingly many results containing the number n, 
eigenvalue of the number operator, have their g— analog with n replaced by [n]. For instance, the 
Hamiltonian of the g— boson oscillator is given by 

^ = E (^^fc + ^"^^fc^fc) = \^^kalak (8) 

k 
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and takes the form 

H = Y,[N] huk (9) 

k 

which has the eigenvalue 

E = ^[n]hujk ^^nhujk. (10) 

k k 

With the inclusion of the zero-point energy, the Hamiltonian, Eq.(|^) will be replaced by the following 
expression 

H = T. I^^fc (2«fe4 + H^k) = E l^^k {[N + 1] + [N]) . (11) 

We note that [A^] is not the number operator. As a consequence, obtaining the distribution function 
in a closed form for the case of the g— deformed oscillators is no longer as straightforward as in the 
case of ordinary boson oscillators. Herein lies the mathematical difficulty, and because of the 
occurrence of [n], obtaining a closed form for the distribution function does not seem to be feasible. 

To deal with this problem there are two methods known in the literature. One consists of 
developing a perturbation theory based on the smallness of the deformation parameter q — 1, or 
rather in powers of 7 = Ing. In other words one can calculate the Hamiltonian, the partition 
function and other quantities of interest valid to some order in 7. This is the method followed by 
Neskovic and Urosevic who derive such perturbation theory results for various thermodynamic 
quantities, demonstrate how to apply this method to some g— deformed boson systems, and study 
typical consequences of small deformations. 

The second method consists of an iteration scheme developed by Song et al which we shall 
use in this work. This method may be briefly reviewed as follows. As a first iteration one may 
approximate the Partition function occurring in Eq.(^) by the zeroth order result, namely the one 
corresponding to the Partition function of undeformed oscillators, or equivalently the zeroth order 
expression for the occupational probability P„. We shall henceforword omit the subscript k. Thus 
we obtain the first order iteration result for the distribution function of the g— deformed harmonic 
oscillator: 

/. = E N Pn = E -T^Pi'' = jz^. E (^"e-'^^ - (12) 
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where x = f3hu. Evaluating the sum and simplifying, we obtain 

^ (g-ie^-l)(ge^'-l) " {qe^^^ - l){q-^eP^^ - I)' ^^^^ 
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This is the form derived by Song et al [Q. The Grand Partition Function which gives rise to this 
distribution function has been studied extensively P, ^ to derive the thermodynamic properties 
of g— bosons as well as in a study of g— phonons. This approximate distribution has the correct 
undeformed limit: 

lim fa = -57 (14) 

It is thus correct in this sense but unfortunately this distribution function does not exhibit the 
correct high temperature limit. We can suspect that this is due to the fact that the above theory 
has discarded the zero point energy. 

Let us first review the well-known high temperature limit of the standard boson oscillators in 
the theory that includes the zero-point energy, where the canonical Partition function is given by 

Z = y e-("+^)- = — = 1-— (15) 

Y 1 - e-^ 2 sinh \x ^ ' 

and the occupational probability is given by 

P„ = = 2e-("+5)a; sinh -. (16) 

Computing the distribution function, we then find 

°° 1 r /I \ 

/ = + 2 = - (- 5: + E ne—) (17) 

and hence 

1 1 _ 1 /e^ + 1\ 
~ 2 ^ - 1 ~ 2 Ve^' - 1/ ■ 
For N non-interacting boson oscillators, the internal energy is times the mean energy of one 
oscillator, thus 

U = Nhuf = Nhu + -^J-^) . (19) 
In the high temperature limit, f3huj << 1, we then obtain 



2 I + phuj + li^huof + 1 

1 1 

N-huj+—{l--(3huo+---)=NkT, (20) 
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the well-known result: the internal energy of quantum boson oscillators is the same as the classical 
value in the high temperature limit. This happens |10] precisely because the term —^NTiuj cancels 
the term arising from the zero-point energy. 



It is evidently necessary to include the zero point energy in the theory of g— deformed boson 
oscillators and we shall proceed accordingly as follows. We see that the distribution function given 
by Eq.(ll) gets modified as 

OO -| 

/. = Eo(h + l] + M)^n, (21) 
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when we incorporate the zero-point energy analogous to the Hamiltonian 

H = hiu{\N + l] + [N]). (22) 

of Eq.(^). These have the correct g — > 1 limits. The occupation probability given by the lowest 
iteration, Eq.(0), corresponding to the undeformed oscillator is thus 

Pi°) = 2sinh-e-("+5)-. (23) 
The basis number can be expressed in terms of g = e'*': 

[n] = ^" - ^7 = (24) 
q — q sinh7 

We can thus express the distribution function as 



^ sinh 77-7 sinh(n + 1)7 . , x _^„,l^^ 

/<7 = E • , ^smh-e ^"+2)-. 
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Evaluating the sum and simplifying the algebra, we obtain 



which can be conveniently put in the form 



(25) 



2{q — q~^) \e'-^ — q — q'^ ~^ — q — q^^ ) ^ 



1 (p2f3hu; _ 1 \ 

f = tl (27) 



4 



or alternatively, 

fi = 2(ge/3^- _ l)(g-ie/5^- - !)■ ^'^^^ 

It is easily verified that in the limit q ^ 1, the above goes over to the form for the standard boson 
oscillator with the inclusion of the zero point energy, given by Eq.(]T8D. It is convenient to rewrite 
the distribution function in terms of partial fractions, thus 

f U I I I ^ (29) 

■''^ 2\l-qe^=' l-g-ig-^ g-^e^ - 1 ge^ - 1 / ' ^ ' 

where Ci + C2 = 1, and 

C^ = ^—^, C2 = (30) 

q — q Q ~ Q 

It is expedient at this point to introduce the chemical potential /i primarily to facilitate some 
manipulations, and we shall set /i = later at the end. Thus we make the replacement 

where z = e^^ is the fugacity. Hence we obtain 

f = -( + ^2 + + 1 f32l 

^ 2 I 1 - qze'P^ 1 - ^e-PE 1 - ^e-^^ ^eP^ - 1 ' ' ^ ' 

We may now introduce the first iteration Partition Function Z which corresponds to this distri- 
bution function. We accordingly determine 

^ = n (1 - ^^'"Y"' (1 - 1^'")"'' - 1)"'"' {^-^-'' - 1)""' (33) 

as the first iteration Partition Function. We may compute the thermodynamic potential for the 
g— boson oscillator as 

VL = --InZ 



- E^{^.ln(l-le-)+C.ln 
+ C2 In (g^e-^^ - l) + Ci In |^ 



1 - ^e'^^ 



^e-^- - 1 ] I . (34) 
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If we now compute the distribution function 



/, = (35) 



we readily reproduce the distribution in Eq. (|32D . At this point we may, for simphcity, return to 
the case where the chemical potential is zero. Thus dropping fi, we obtain the distribution function 



fg- 2^.(^1 



Ci , C2 . C2 Ci \ 

+ 73^?^ T (36) 



- qe-P^ 1 - g-ie-/3^ q-^e^^ - 1 qe(^^ - 1 



and the Partition Function 



Z = n (1 - q-'e^^Y'"'' (1 - ge^^)"^""' (ge"^^ - l)"^""' (g-^e"^^ - l)"^""' . (37) 

E 

It is readily verified that these forms possess the correct g — > 1 limits, Eq.(|18]) and 

^=n^. (38) 



We may now compute the internal energy: 



f/ = -|ln2 (39) 

= ^ V p f I C2 Ci \ 

2^ \qe-f^^ - 1 q-^e-f^^ - 1 1 - q-^ef^^ 1 - qe^E J ' ^ ' 

which clearly has the correct g — 1 limit, namely 

(41) 



It is convenient to express the internal energy in an alternate form, since the internal energy is 
times the mean energy of one oscillator, 

1 / C C C C \ 

U = -NEi ^ + + , ' ^ + ' ^ , (42) 



The undeformed limit is 



U,_, = J-NE ^^^^ = iiVsfl+i (43) 

^"^ 2 l-cosh/?E 2 e^'-l ^ ^ 
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as expected. For g 7^ 1, we may cast it in the form 

U = -Ne'^^ ( ] . (44) 

g-g-i Vl + g2-2gcosh/3E 1 + q-^ - 2q-^ cosh f3 E J ^ ' 

Upon introducing q = e"' , we derive the result in a convenient form 

U = NE = '-NE ^^^hf coshf 

2(cosh/3i?-cosh7) 2 sinh^sinhto ^ ^ 

which has the famihar g ^ 1 hmit, namely 



Let us first consider the low temperature limit. For PE >> 1, the internal energy will not have 
the correct low temperature limit for arbitrary q since it will depend on 7. We shall thus consider 
small deformations only, 7 << 1, or rather 7 << jSE in terms of the dimensionless quantity. For 
very low temperatures, we thus have l3E >> 1 >> 7. In this case then 

1 ,.^sinh/5E/2cosh/5E/2 1,^^ 

Ug^oo = m -NE , ; ^ , ' = -NE 47 

/3-.00 2 (sinh/5E/2)2 2 ^ ' 

which is the zero-point-energy. The low temperature limit is the pure quantum effect as expected. 
Thus the theory at T = is the same as the undeformed theory W^ . 



Next we examine the high temperature limit. In the case g = 1 , it has the correct classical 
limit, namely 

Uq=i NkT, (3E«l (48) 

but this is not the case when q ^ 1 for arbitrarily large values of 7 7^ 0. We accordingly consider 
small deformations defined by 7 << (3E. At high temperatures we thus have 7 << (3E « 1. 
Consequently we find from Eq.(H^) that 



Ug^, = \nE^ = NP, (49) 

which is the expected classical limit. This implies that as long as the deformations are small, the 
high temperature limit of the g— deformed boson oscillators is the same as that of the undeformed 
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theory. It is gratifying to note that this is parallel to the expectation that at T = 0, the theory of 
the g— deformed oscillators is the same as that of the undeformed theory |[TT|] . 



In summary, we have developed the theory of the statistical mechanics of g— deformed boson 
oscillators incorporating the zero-point-energy. We have confronted the theory against the correct 
limits at very high temperatures and very low temperatures, we thus conclude that (a) At very low 
temperatures, T = 0, the theory is no different from the undeformed theory of boson oscillators and 
(b) at very high temperatures, T = oo, the internal energy corresponds to the classical limit and 
the theory reduces to the undeformed case at T = oo and (c) these correct limits are possible only 
after incorporating the zero-point energy in the spectrum. 

The theory of g— deformed fermion oscillators pose different challenges. This investigation will 
be reported in a future publication. 

Acknowledgments 

I am grateful to P. Quarati and A. Lavagno of Politecnico di Torino, Torino, Italy for the kind 
hospitality where this work was begun. I am thankful to A. Lavagno for valuable discussions and 
drawing my attention to some papers on the subject of g— oscillators. I thank R. Ramachandran, 
director of the Institute of Mathematical Sciences (IMSc), Chennai, India, for the hospitality ex- 
tended to me, where I did much of this work during part of my sabbatical leave in 1997. I would 
like to thank colleagues at IMSc: G. Rajasekaran for many valuable discussions and for educating 
me on the nuances of the g— deformation algebra; R. Parthasarathy for discussions on the general 
subject of g— oscillators, and R. Jagannathan and G. Baskaran for stimulating discussions. 

REFERENCES AND FOOTNOTES 



1. L. Biedenharn, J. Phys. A22, L873 (1989) 

2. Q.Yang and B.Xu, J. Phys. A 26, L365 (1993); N. Aizawa, ibid A26, 1115 (1993); K. 
Viswanathan et al ibid, A25, L335 (1993); E. G. Floratos, J. Phys. A24, 4739-4750 (1991); 
T. Altherr and T. Grandou, Nuc. Phys. B402, 195-216 (1993) 

3. R. Parthasarathy and K. S. Viswanathan, J. Phys. A 24, 613-617 (1991), K. S. Viswanathan, 
R. Parthasarathy and R. Jagannathan, ibidA25, L335-339 (1992). 



8 



4. A.K. Mishra and G. Rajasckaran, Mod. Phys. Lett. A 9, 419-426 (1994); A. K. Mishra and 
G. Rajasekaran, Pramana 45, 91-139 (1995). 

5. H.Song, S.Ding and I.An, J. Phys. A 26, 5197 (1993). 

6. G. Kaniadakis and P. Quarati, Phys. Rev. E49, 5103 (1994); G. Kaniadakis, Phys.Rev. E49, 
5111 (1994); G. Kaniadakis, A. Lavagno and P. Quarati, NucL Phys. B 466 527 (1996). ibid, 
Phys. Lett. A 227 227 (1997). 

7. P. Neskovic and B. Urosevic, Int. Journal Mod. Phys.AAl, 3379 (1992) 

8. P. Narayana Swamy, Int. Journ. Mod. Phys. B 10, 683-699 (1996). 

9. P. Narayana Swamy, Mod. Phys. Lett. B 10, 23-28 (1996). 

10. W. Greiner, L. Neise and H. Stocker, Thermodynamics and Statistical Mechanics, Springer- 
Verlag Inc. (New York) 1995. 

11. J. L. Birman, Phys. Lett. A 167, 363 (1992). That deformation does not prevail at T = 
is also apparent in a g'— deformed theory of fermionic oscillators. See e.g., G. R. Lee and J. P. 
Yu, Phys. Lett. A 164, 164 (1992). 



9 



